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Abstract 

We consider constraints from flavor physics on two-Higgs-doublet models (2HDM) with 
general, flavor-diagonal, Yukawa couplings. Analyzing the charged Higgs contribution to 
different observables, we find that b — ► S7 transitions and AMs d restrict the coupling 
Xu of the top quark (corresponding to cot/3 in models with a Z2 symmetry) to |Aff| < 1 
for rrifj+ ^ 500 GeV. Stringent constraints from B meson decays are obtained also on 
the other third generation couplings A&& and A rr , but with stronger dependence on m H +. 
For the second generation, we obtain constraints on combinations of X ss , A cc , and A^u 
from leptonic K and D s decays. The limits on the general couplings are translated to 
the common 2HDM types I - IV with a Z2 symmetry, and presented on the (m#+,tan/3) 
plane. The flavor constraints are most excluding in the type II model which lacks a 
decoupling limit in tan/3. We obtain a lower limit > 300 GeV in models of type II 

and III, while no lower bound on is found for types I and IV. 
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1 Introduction 



The LHC will hopefully shed some light soon on the mysteries of electroweak symmetry break- 
ing. Meanwhile we are free to speculate whether nature is endowed with new physics beyond 
the minimal Higgs framework with one scalar doublet, as incorporated in the standard model 
(SM). One of the simplest - and probably the most thoroughly studied - extension of the Higgs 
sector is the two-Higgs-doublet model (2HDM), in which the SM Higgs sector is duplicated 
by introduction of another doublet of SU(2)l [1]. The primary motivation for considering the 
2HDM is often sought in supersymmetry [2], since the minimal supersymmetric standard model 
(MSSM) requires at least two Higgs doublets in order to cancel gauge anomalies and keep the 
superpotential analytic. This second requirement fixes the pattern of Yukawa couplings at 
tree-level. Many interesting results [3-8] have also been obtained for more general two-Higgs- 
doublet models. These models can be thought of as effective theories for supersymmetric Higgs 
sectors beyond the MSSM (see e.g. [9] and references therein). In addition they provide a 
generic description for UV extensions of the SM that do not involve supersymmetry, but still 
contain two Higgs doublets in the low-energy effective theory. Two-Higgs-doublet models of 
the general form are the focus of the present work. 

The additional scalar doublet in the 2HDM introduces new and interesting phenomenology, 
such as new sources of CP- violation, flavor-changing neutral currents (FCNC), the presence of 
a charged Higgs boson, and possibly even a dark matter candidate [10]. In the Yukawa sector, 
the additional doublet leads to a large number of couplings between fermions and scalars which 
are not fixed by the mass generation mechanism. When FCNCs are absent at tree-level, the 
most striking phenomenological difference between the SM and the 2HDM Yukawa sectors is 
the presence of a charged Higgs boson and its associated charged current. It is mainly through 
this interaction the 2HDM can reveal its presence in low-energy observables, such as the decay 
rates of K, D, and B mesons, complementing searches for charged Higgs bosons at high-energy 
colliders. Results on the 2HDM Yukawa sector are therefore useful in both these contexts. 

Constraints on the 2HDM parameters have been extensively studied in the literature, in 
particular for the type II model [11] which is the form of the Yukawa sector in the MSSM at 
tree-level. For recent studies discussing flavor constraints on the charged Higgs boson in the 
MSSM, see [12]. To differentiate between the MSSM Higgs sector and a more general theory 
with two Higgs doublets, it is interesting to identify differences in phenomenologically viable 
signals, such as a charged Higgs boson contribution to flavor physics observables. Combining 
future LHC results with flavor physics could also be a means of testing the Yukawa couplings 
to more than one type and generation of fermions. To determine the reach in the flavor sector, 
we extend the previous results for the 2HDM type II to flavor constraints in the more general 
setting of arbitrary flavor-diagonal Yukawa couplings, i.e. with no assumptions of coupling 
universality between fermions of different types or generations. Limits on the relevant couplings 
are evaluated for each observable separately. As a special case, we then translate the general 
results into constraints on the (m#+, tan/3) plane for types I - IV [13] of the 2HDM with a Z 2 
symmetry^] 

1 The phenomenology of the more unusual models type III and IV, which have different couplings for the 
down-type quarks and the leptons, has recently been studied in [14]. 
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The organization of this paper is as follows: in Section [2] we give a brief introduction to the 
2HDM, describing how the parameters of the potential are constrained by theoretical arguments 
and electroweak precision data. Following this discussion of the scalar sector, Section [3] deals 
with the 2HDM Yukawa sector in some detail. The different flavor observables of interest are 
introduced in Section HJ where also also our main results are presented as constraints on the 
individual Yukawa couplings from each observable separately. In Section [5] we specialize the 
general results of Section H] to models with a Yukawa Z 2 symmetry. Finally, Section [6] presents 
the conclusions of this work. 



2 2HDM Fundamentals 

In a generic basis the most general, renormalizable, Higgs potential for two identical doublets 
$i and $2, with hypercharge Y — 1, is given by [4] 



^2HDM = m\ x <b\&\ + $ 2 $2 ~ m{ 2 §[®2 + h.C. 

+ (V^l) 2 + ^A 2 ($ 2 $ 2 ) 2 + A 3 (V^l) ($ 2 $ 2 ) + A 4 (V^) ($2$! 

+ |^A 5 (V^)' + [A 6 (V^x) + A 7 ($J$ 2 )] ($t$ 2 ) + h.c. 

In total there are 14 free parameters in this model, including the complex phases which may be 
present in A5 6 7 and m\ 2 . Restricting ourselves to the case without CP violation, all parameters 
are assumed to be real, and the number of free parameters is reduced to 10. For the electroweak 
(EW) symmetry to be broken, the scalar mass matrix has at least one negative eigenvalue. 
At the minimum, m\ x and m 22 can be eliminated in favor of the vacuum expectation values 
= Vi/y2. The 2HDM is invariant under a unitary transformation of the scalar doublets 
[5,6]. A basis in the doublet space is chosen by specifying tan/3 = v 2 jv\. The overall scale (like 



2 <ist, 



, 2 <f>t, 



in the SM) is set by v z 



(246 GeV) . For the sake of simplicity, we take A6 = A7 = 



in this section. This is required if the potential should respect a discrete Z 2 symmetry imposed 
on the $i fields^] We will return to the importance of discrete symmetries later. The number 
of free parameters, v being fixed, is down to six (plus tan/3 which specifies the basis). 

When the EW symmetry is broken, three of the eight scalar degrees of freedom are used for 
the masses of the gauge bosons. The physical spectrum thus contains three neutral Higgs bosons 
and the charged Higgs pair H ± . Two of the neutrals are CP-even (h and H, with rrih < ma), 
and one (A) is CP-odd. The CP-even scalars mix with an angle a. At the minimum, the 
original doublets are expanded according to 



1 / V2 (G + cos /3 - H + sin (3) 

v 2 \ v C os (3 — h sin a + H cos a + i (G° cos (3 — Asm (3) 



(2) 



2 The Zi symmetry can still be (softly) broken by dimension two terms for non-zero values of m\ 2 - 
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1 / V2 (G + sin (3 + H + cos (3) 

$2 = -== V 7 | • (3) 

v 2 I s in (3 + /i cos a + H sin a + i (G° sin /3 + A cos /3) 

Using the notation sp = sin /3, eg = cos /3, the tree-level Higgs masses are given in terms of the 
remaining potential parameters of Eq. (CQ) by 



m\ = ^- \,v 2 (4) 



2 



for the CP-odd Higgs, 



m 



H,h 



Ml + M 2 2 2 ± J {Ml - M| 2 ) 2 + 4 (M 



2 

for the CP-even Higgses, where the mass matrix is 



f 2 ) 2 



(5) 



and finally 



-spcp £p J \ (A 3 + A 4 )s /3 c /3 A 2 s^ + \ 5 c 2 

m 2 H+ = m\ + t^ 2 (A 5 - A 4 ) (7) 

for the charged Higgs mass. Working at tree-level, the mass relations allow the input parameters 
Ai_5 to be substituted by the four Higgs masses and the neutral sector mixing sin(/3 — a), giving 
a set of physical input parameters. 

The parameters in Eq. (TjQ) are bounded by the requirement that the potential is stable 
[15], i.e. that no direction in field space exists for which V — > —oo. Necessary and sufficient 
conditions for stability are given by 



Ai > 0, A 2 > 0, A 3 > — \AiA 2 , 

A 3 + A 4 — |A 5 | > — a/AiA 2 . 

The parameters can also be bounded from above by arguments of perturbativity and pertur- 
bative unitarity of the 5-matrix for Higgs and longitudinal vector boson scattering [16]. We 
consider the upper limits |Aj| < 47r (perturbativity) and \L\ < I6n for the eigenvalues L of the 
S'-matrix. Assuming the 2HDM is perturbatively defined up to some scale A, these conditions 
should in principle be fulfilled at all scales fi < A. The resulting limits on the A« can then be 
studied in an RGE-improved treatment [17]. Since this is not the primary aim of this paper, we 
limit our discussion of the theoretical consistency to the input (EW) scale. When the potential 
has an exact Z 2 symmetry (m\ 2 = 0); the requirement of stability, perturbativity, and unitarity 
leads to upper limits on the Higgs masses of m h < 600 GeV, m H < 870 GeV, m A < 870 GeV, 
and rriH+ ^ 780 GeV. Breaking the Z 2 symmetry softly, the upper limits are further increased, 
as should be expected when introducing the new scale m\ 2 . In the decoupling limit [4], the 
Higgs masses scale like rrih ~ v, and m H A,H+ ~ m s (1 + v2 / m2 s), where ms is some high scale 
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Figure 1: Mass differences for 2HDM points which give contributions to S, T, and U in agree- 
ment with the present limits [20] at the 2 a level. 

ms > w at which the heavy doublet is integrated out. The existence of the decoupling limit 
shows the absence of an absolute upper limit for the heavy Higgs masses in the general 2HDM. 
We should therefore focus on deriving lower limits on the masses, which will require additional 
input from experiments. 

Before considering the Yukawa sector, we want to briefly discuss electroweak precision tests 
and the 2HDM contribution to the oblique parameters S, T, U [18]. Fixing rrih = 114 GeV, 
we perform a scan over the other 2HDM masses (with > rrih) and sin(/3 — a). For each 
point, the results of [19] are used to evaluate the oblique parameters. Points which fulfill the 
PDG limits [20] at the 2 a level are illustrated in Fig. [TJ The figure shows the mass splittings 
\rriH+ — tua\ and \uih+ — Tns\, where m| = m 2 H sin 2 (/9 — a) + m\ cos 2 (/5 — a) is a combined scalar 
mass. We see that both these mass splittings cannot be simultaneously 0(v ) for the 2HDM 
to be compatible with the EW precision tests. The limits m H + ~ (independent of (3 — a), 
rriH+ — ran (with cos(/3 — a) =0), and m H + ~ rrih (with sin( / 9 — a) = 0) correspond to the 
known cases when custodial symmetry ensures T ~ [8,21]. 

3 The Yukawa sector 

Assuming CP conservation, the general Yukawa Lagrangian of the 2HDM, expressed in the 
fermion mass eigenstates, is [6] 

-CY^LQivYUR + Q&rtfDR + Q&tfLR + Q^ (8) 

where $j = io"2$*. The rjf, with F = {U,D,L}, are real 3x3 Yukawa matrices, subject to 
the requirement that masses for the fermions are generated by the combination coupling to the 
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vevs: 



M 1 



V2 



(rj F cos j3 + rj F sin /?) 



(9) 



M F is here the (real and positive) mass matrix. It is convenient to introduce the short-hand 
notation 

k f = r) F cos (3 + rj F sin (3. (10) 

The combination orthogonal to n F , and which does not participate in giving mass to the 
fermions, is then 

p F = —Tj F sin (3 + r/f cos (3. (11) 

Performing a change of basis on the Higgs fields, we can also re-express the Yukawa Lagrangian, 
Eq. ([H]), in the basis (Hi,H 2 ) where only one Higgs doublet (here taken to be Hi) acquires a 
vev. The resulting Lagrangian becomes 



-Cy = Q L HiK Ur+Q l Hik Dr+Q l Hik Lji+Q L H 2 p Uji+Q L H 2 p D^+Q L H 2 p Lji^ (12) 

H\ (with vev V\ = v) generates the fermion masses and consequently couples with n F , while 
if 2 has v 2 = and couples with p F . 

The Yukawa matrices for the Higgs doublets determine the couplings of the physical scalars 
to the fermions, which in an arbitrary basis are given by 



s/2 



K U Sp- 



+ P Cf3_ a 



Dh H —D 

Uh + -^=U 
y/2 



P U s P - 



DH+—D l5 p D DA 
UH - ^=U l5 p u UA 



V2 



K L S(S- a 



Lh + -^-=L 

V2 



(13) 



K L C/3- 



p L sp- a 



+ 



+ P C(3- a 

U{V CKM p D P R - P u V C kmPl)DH + + Vp L P R LH + + h.c. 



LH + -^L^p L LA 
V2 



The general form of Eq. (1T31) exhibits tree-level FCNCs unless the p F are diagonal. A sufficient 
condition for this to occur is that each fermion type F = {D, U, L} couples to only one Higgs 
doublet [22]. In Eq. ([8]) this condition translates into demanding either t][ = or 7] 2 = 0, which 
through Eqs. ffTDl . ffTTl) leads to the relations p F = k f cot [3 and p F = — K F tan/5, respectively. 
A way to accomplish vanishing rjf naturally is to impose a symmetry to prevent some of the 
couplings from appearing in Eq. , for example a Z 2 symmetry under which one Higgs doublet 
and some of the right-handed fermion fields are odd. Models with such a symmetry are called 
2HDM types. We choose a convention for the types I - IV given by [13]jf| The assignment of 
Z 2 charges, and the emergent coupling relations, are specified in Table [H 



3 The notation and meaning of the different types varies in the literature. Sometimes type Y (III) and type X 
(IV) is used. In supcrsymmctry, type III usually refers to the general model encountered when the Zi symmetry 
of the tree-level type II model is broken by higher order corrections. 
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Type 


Ur 


Dr 


Lr 


p u 


P D 


P L 


I 


+ 


+ 


+ 


k u cot (3 


K D cot (3 


K L COt /? 


II 


+ 






k u cot (3 


—k d tan/3 


— k l tan /3 


iii/y 


+ 




+ 


K U cot (3 


— k d tan/3 


k l cot (3 


IV/X 


+ 


+ 




n u cot (3 


k d cot (3 


—k l tan /3 



Table 1: Assignment of Z 2 charges for the right-handed fermions, and the resulting relations 
among Yukawa coupling matrices in the Z 2 -symmetric types of 2HDM Yukawa sectors. The 
Higgs doublets have Z 2 quantum numbers — ($1) and + ($2)- 



The ^-symmetric 2HDM types constitute minimal models for the Yukawa sector in the 
sense that they are one-parameter descriptions. From a phenomeno logical viewpoint, a rather 
natural step beyond the types are models with Yukawa alignment [23]. Assuming that the 
Yukawa matrices for each type of fermions are proportional by a constant £ F , one can define 
three tan/3-like parameters through p F = tan (3 f k f , where tan/3 F = ■ Since k f and 

p F can be diagonalized simultaneously, the model is free of tree-level FCNCs. A further gen- 
eralization is offered by models where the off-diagonal Yukawa elements are non-zero, but still 
sufficiently small to avoid the experimental bounds from FCNC Using a a Hermitian six-zero 
texture for the Yukawa matrices ll a phenomenologically viable ansatz for this type of models 
was introduced by Cheng and Sher [25] 



[p F h = t ] v ~ (™ sum). (14) 

It leads to a natural suppression of the off-diagonal elements for A ~ 0(1) by the quark mass 
hierarchy, in particular for the first generation where the experimental constraints are most 
restrictive. Different scenarios have been investigated in [26]. The Cheng-Sher ansatz provides 
a parametrization useful also for general, flavor-diagonal, Yukawa couplings. From the definition 
of k, f , it follows that the diagonal elements are 



72 



rrii 



[p F }u = [\*']u V "J"* = [X F }ii[n F }u (no sum). (15) 

This relation describes general flavor-diagonal couplings with a preserved mass hierarchyjf] The 
nine (one for each fermion) are free tan/3-like parameters. Both the case of alignment (Af x = 
X22 = A^ 3 = tan/3^ 7 , similarly for D and L), and the models with a Z 2 symmetry (e.g. type II: 
\V = cot /3, A^ = X F = — tan (3) are contained as special cases of this general framework. We 
shall use this phenomenological approach in the following. The Yukawa couplings have been 
implemented in the generic code 2HDMC [27] used throughout this work. When there is no 



4 The phenomenology of the related case with four-zero textures was recently discussed in [24]. 
5 From the appearance of equation (|15|) one could refer to this case as generalized alignment, although the 
motivation in terms of a proportionality between the Yukawa matrices no longer holds. 
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risk of confusion, we use the notation X a = A33, X bb = A33, etc. Inspired by the ^-symmetric 
models, we further restrict the couplings to, at most, allow for hierarchies of the order 

m d \X dd \ ~ m u \X uu \ ~ 

m s |A ss | ~ m c |A cc |, m c \X cc \ ~ m b \\ bb \, m b \X bb \ ~ m t \ X tt \ (16) 
m e \X ee \ < m^lA^I < m r |A TT |. 

This means couplings to the first generation are always neglected. Couplings to the nearest 
neighbors in mass can have a hierarchy allowing them to be of similar magnitude, hence they 
must be considered simultaneously in observables where both couplings appear. No strong 
hierarchy is allowed between generations of equal charge fermions, since this would be unnatural 
under the Cheng-Sher assumption. 

4 Constraints from flavor physics 

The most stringent experimental constraints on the 2HDM come from flavor physics. When 
the Yukawa couplings are diagonal in flavor space, one could still have sizeable contributions 
to many observables due to charged Higgs exchange, either at tree-level or at one-loop. We 
shall focus on such processes here. The Yukawa coupling dependence from the lowest order 
charged Higgs contribution to the interesting observables is listed in Table [2J together with the 
current experimental results and the SM expectations. For the numerical evaluation we use 
Super Iso v2.6 [28], both for the SM and the 2HDM contributions. The slight differences in some 
of the SM values compared to earlier published results are explained by parametric updates. 
The top quark pole mass is set to the latest combined Tevatron value m t = 173.1±0.6±1.1 [29], 
corresponding to m t = 164.0 GeV using the 2- loop relation. For the other quark masses, we 
use the PDG MS values [20] m b = 4.20^;^ GeV, m c = 1.27l^ GeV, and m s = 104^ MeV 
as input. When calculating leptonic- and semi-leptonic decays of mesons, the running quark 
masses appearing in the H + couplings are evaluated at the scale of the meson mass. This 
ensures that the effective four-fermion operators mi^^r^l, obtained when integrating out H + , 
renormalize in the same way as the corresponding quark currents [30]. Constraints on the 
2HDM parameters from each observable are presented as exclusion regions at 95% C.L., using 
Gaussian two-sided confidence intervals and one degree of freedom. 

Negative searches for the charged Higgs at LEP lead to a lower limit m H + > 78 - 90 
GeV [36]. The exact value depends on the H + decay mode. Strictly speaking, these limits 
are only valid under the assumption that the H + decays with total BR = 1 in the modes 
H + — > t + u t and H + — > cs, with certain benchmarks investigated for H + — > A showing 

similar exclusion. Even if the 2HDM with general Yukawa couplings could certainly allow for 
exotic decay channels, hierarchical couplings still favor decays into these channels. We will 
therefore respect the LEP limits in the following and not consider the possibility < my/. 

The first flavor observable we consider is the branching ratio for the rare FCNC inclusive 
decay B —>■ V s 7, which receives contributions from the charged Higgs at the same level as the 
W contribution in the SM. Using an effective Hamiltonian approach, the leading order 2HDM 
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Observable Couplings Experimental value Standard Model 



BR(B - 


y X s 7) 


tt 5 *™ 


(3.52 ± 0.23 ± 0.09) x 10~ 4 [31] 


(3.07 ± 0.22) 


x 10~ 4 


A (B -> 


• K*i) 


\-H-\hh 


(3.1 ± 2.3) x 10~ 2 [28] 


(7.8 ± 1.7) x 


10~ 2 


AM R 




A 2 , 


(0.507 ± 0.004) ps" 1 [31] 


(0.53 ± 0.08) 




BR(B U - 


-> rz/ r ) 


Xbb^TT 


(1.73 ±0.35) x 10~ 4 [32] 


(0.95 ±0.27) 


x 10~ 4 


£dIv 




Abf,A Tr , A cc A rr 


0.416 ±0.117 ±0.052 [33] 


0.30 ±0.02 




R m (K- 


-»■ /^V) 


Ass Xfxfj, 


1.004 ±0.007 [34] 


1 




BR(D S - 




XssXfxfxi A CC A^ 


(5.8 ± 0.4) x 10" 3 [35] 


(4.98 ±0.15) 


x 10~ 3 


BR(D S - 


-> rz/ T ) 


AssA rr , A cc A rr 


(5.7 ±0.4) x 10~ 2 [35] 


(4.82 ±0.14) 


x 10~ 2 



Table 2: List of observables sensitive to a charged Higgs contribution at lowest order, their 
dependence on the 2HDM Yukawa couplings, the current experimental and SM values. For 
cases where the experimental value is a combination of several values, we reference the paper 
where the combination is performed. The SM values are evaluated with Superlso v2.6 [28], 
using the methods and input parameters described in the text. 

contribution to the relevant Wilson coefficients at the matching scale fiw is [37] 

5CVT M (M = l^(y) - X tt X bb F^(y) (17) 

where y = m 2 / m 2 H+ , and 

7 {y) ~ 240,-1)3 + 4(y-l) 4 [Uy 

8 {y) ~ 8(y-l)3 4(,-l) 4lny 
U2) ( » y(3-5y) y(3y-2) 

The lowest order expression given by Eq. (1171) captures the essential dependence on the Yukawa 
couplings A^ even when higher order corrections are included. The numerical evaluation 
of BR(B -> X s -f) [28] is performed to NNLO for the SM [38], and to NLO for the 2HDM 
contribution [39]. Results on the excluded regions of parameter space (X tt , Abb) are shown in 
Fig. [2] for two values of the charged Higgs mass: m#+ = 150 GeV and m#+ = 400 GeV. These 
were chosen having the well-known limit m H + > 300 GeV for the 2HDM II in mind. When 
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Figure 2: Constraints on (X tt , X bb ) from BR(B — > X s 7) for fixed = 150 GeV (left) and 

m H + = 400 GeV (right). The highlighted region is excluded at 95% C.L. 

m,H+ = 150 GeV, we find that a large fraction of the parameter space is excluded, with a 
preference for smaller values of the couplings in the allowed region. We also see allowed regions 
in Fig. [2] with simultaneously large and equal sign values for A&& and \ tt - For these fine-tuned 
values, a cancellation occurs in the new physics contribution to the Wilson coefficients. The 
situation is similar for m H + = 400 GeV, except that the size of the allowed region increases. 
As expected from Eq. (TPTI) . taking A&& — > is in general not enough to avoid generating a large 
H + contribution, but Atf — > always is. 

Another quantity measured in b — > 57 transitions is the degree of isospin asymmetry in the 
exclusive decay mode B — > K*y, defined as [40] 

T(B° -> if *° 7 ) - r(g- -> A-*- 7 ) , . 

°- r(B°^i?*°7) + r(B-^i?*-7)' 1 ' 

This observable gives more stringent constraints on the model parameters than BK(B — > X s y) 
in several MSSM scenarios [41]. Using the NLO prediction we investigate, for the first time, 
constraints from A _ on the Yukawa sector of the general 2HDM. The excluded region in 
(A«, At;,) is presented in Fig. [3l To facilitate a comparison with the results for the branching 
ratio, the same values as in Fig. [2] are chosen for m#+ . From Figs. [2] and El one notices a 
similarity in the regions excluded by BK(B — * X s y) and Ao_. This results from the common 
dependence of both observables on the Wilson coefficient CV: BK(B —* X s y) is proportional 
to Cj, while Ao_ varies like Cj 1 to first order. The isospin asymmetry results exclude the 
large-coupling solution observed in Figure [2] for same-sign couplings. 

The AMs d and AM# s mass differences measured in B°-B° mixing are sensitive to charged 
Higgs exchange through box diagrams involving internal top quarks. Including the leading 
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Figure 3: Constraints on (X tt , A^) from the isospin asymmetry A _ for fixed m#+ = 150 GeV 
(left) and m H + = 400 GeV (right). The highlighted region is excluded at 95% C.L. 



order contributions, the 2HDM prediction for AM# d is [42,43] 



AM, 



G' F mif^B d M B \V t * d V tb \% 
24vr 2 



lww(y w ) + l WH (y w , y 11 , x) + i HH (y H ] 



(20) 



where y l = ml/mf {i = W,H + ), x = m 2 H+ /m w , and 



6 



i ww 



Iwh — x 2 tt y H 



(i - y w ) 2 

(2x-8)lny^ 

{l-x){l-y H f ^ (i_ x )(i_ y w)2 - {l _ y w ){l _ y H ) 



6 



y 



w 



w 



\i-y 

6x\ny w 



hi y 



w 



-2y 



(21) 



4 „.H 



K t y 



H\2 



+ 



2y In y 

(l-2/ H ) 3 



Approximate effects of short- distance QCD corrections are incorporated in Eq. ( 1201) through 
the factor rjb = 0.552 [44]. The non-perturbative decay constant /# d and the bag parameter B^ 
are evaluated simultaneously from lattice QCD. We use the value fB d B d = 216 ±15 MeV [45], 
adding all errors in quadrature. Since the uncertainty in fB d B d is correlated to the error in 
the corresponding parameters for B s decays, and since the theoretical uncertainty dominates 
the experimental uncertainties, AM^ and AM# S do not provide independent constraints on 
the 2HDM. We therefore consider only AM^ which has the smallest total uncertainty of the 
two. In Fig. H]the excluded region for |A^| is shown as a function of m#+. Nothing can be said 
about the sign of X u since it enters AM^ quadratically. From Fig. H] it is clear that large values 
of | A f t | are disfavored; we find that \X U \ > 1 is ruled out for m H + < 500 GeV, independently of 
the other couplings. 
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Figure 4: Constraints on (m H +, \ \tt\) from AM Bd 
C.L. 



The highlighted region is excluded at 95% 



There are a number of decays which the charged Higgs can mediate already at tree-level. 
The first such observable we consider is B u — > rz/ T , which has a small branching ratio in the SM 
caused by helicity suppression. Since no suppression occurs for the H + , the two contributions 
may be of similar magnitude, leading to sizable interference and a reduced branching fraction. 
The 2HDM contribution factorizes, giving [46] 



BR(i? 1J 



T1S T 



with 



R 



2HDM 
Btu 



8vrr B 



1 - 



-m B m 



mz 



B 



2HDM 



m 



B 



Til 



Aw,A T 



(22) 



(23) 



For the B decay constant we use the value fs = 190 ± 13 MeV [45], which together with V u b 
constitutes the main theoretical uncertainty. We use the combined CKMfitter value \V u b\ = 
3.87 ± 0.09 ± 0.46 [32]. Since only the product XbbKr enters Eq. (123]) . the exclusion region is 
presented for this quantity versus m H +. The result is shown in Fig. [5], where \ bb \ TT /\\ b b\ TT \ 1 / 2 
on the y-axis corresponds to tan/3 in the familiar type II model (for positive values). B u — > tv t 
excludes a large region of the parameter space, with stronger limits obtained on the product 
of couplings for lower values of m#+. Around m#+ = 100-250 GeV, there is a region for 
X b bX TT > 20 where the interference is exactly twice the SM amplitude, but opposite in sign, 
leading to a cancellation in Eq. fl2"31 . 

Another third generation observable, complementary to B u — > tu t , is the semi-leptonic 
mode B — > Dru T . With the 2HDM contribution included, the differential decay rate with 
respect to w = vb • Vd is given by [47,48] 



dT(B -> Dtv t ) 
dw 



V cb \ 2 m\ 



1927T 



-pv{w)x 



mz 
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B 



m b X b b - m c X cc 

H w ) 7 T^ At 

(m b - m c )m A H+ 



ps{w) 



(24) 
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Figure 5: Constraints on (m H +, \ bb \ TT ) from BR(B U — > ru T ). The highlighted region is excluded 
at 95% C.L. The y-axis normalization corresponds to tan/3 in the 2HDM II (for positive values 
only). 



where py {ps) are vector (scalar) form factors and t(w) = m 2 B + m 2 D — Iwm^mB- It has been 
shown [48] that comparing differential distributions directly would be a superior method to 
extract the charged Higgs contribution. However, since the collected statistics is still too low, 
this method has so far not been pursued experimentally. To reduce the uncertainty from the 
vector form factor, we consider the ratio [47] 

_ BR(£? -> Dru T ) 
imv ~ BR{B - Dev e ) ' (25) 

where the 2HDM contributes only to the numerator. The resulting constraints based on ^uiv 
are shown in Fig. [6] for m#+ = 100, 200 GeV. Similarly to the result for B u — > tu t , the 2HDM 
contribution in Eq. (1241) can be twice the SM contribution with opposite sign, leading to the 
two disjoint exclusion regions as observed in Fig. When A cc <C A&& - like in the 2HDM II at 
high tan/5 - the effective constraint on A^A TT can be combined with that from B u — > tv t to 
cover fully the cancellation region observed for low m#+. 

The last B decay we consider is B s — > p + fi~, which has so far not been observed experi- 
mentally. The SM prediction for the branching ratio is 

BR(B S -> p + fi~)sM = (3.2 ± 0.5) x 10~ 9 , (26) 

while the current experimental limit, derived by the CDF collaboration, is [49]: 

BR(B S -> p + fi~) < 5.8 x 10" 8 (27) 

at 95% C.L. The gap between the SM prediction and the current experimental limit makes 
this observable particularly interesting in SUSY, since this difference leaves room for SUSY 
contributions. In the 2HDM however, we found that the experimental limit can be reached only 
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Figure 6: Constraints on (A CC A TT , X hh \ TT ) from BR(£> — > Dtu t ) for fixed m#+ = 100 GeV (left) 
and m H + = 200 GeV (right). The highlighted region is excluded at 95% C.L. The axes are 
normalized such that the values correspond to tan/3 (cot (3) in the Z 2 symmetric versions of the 
2HDM. 



for very large values of A w and A&& or X u and small charged Higgs masses. The constraining 
power of B s — > fi + Lt~ in this study is hence rather limited as compared to the other flavor 
observables and therefore we do not investigate it further. 
Recently it was suggested [34] to use the observable 



R 



123 



V us (K i2 ) Kd(0+ 0- 



V US (K, 



("3 J 



(28) 



from K — > ixv ^ transitions to search for charged Higgs contributions. In this expression V US {K^) 
means V us measured in leptonic kaon decay with i particles in the final state (two leptons and 
i — 2 neutral pions). The same notation applies for 7i> 2 while V u d(0 + — > + ) refers to V u a 
measured in nuclear beta decay. The 2HDM contribution to Rnz is given by 



R 



(23 



K+ 



M 2 m 



(29) 



where we use m^/m s = 1/20 [20]. The major source of uncertainty in extracting i^ 23 exper- 
imentally originates in the necessary input of the form factor ratio fx/fw from lattice QCD. 
From an unquenched calculation with staggered fermions, the HPQCD and UKQCD groups 
have determined fx/fn = 1.189 ±0.007 [50] which is the number adopted in [34]. We will make 
use of that same result, keeping in mind that alternative determinations with larger errors ex- 
ist. Using these would of course reduce the constraining power of K — > iiv^ transitions. Fig. [7] 
shows the result for the product A as A w as a function of m#+. Exclusion is obtained for low 
rriH+ and large absolute values of the coupling, (A^A^I > 20 for m H + = 100 GeV. The i?£ 2 3 
observable is interesting since it provides constraints on the Yukawa couplings to the second 
generation of fermions. 
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Figure 7: Constraints on (m H +, A SS A MM ) from R^2z{K —> 
eluded at 95% CL. The axis normalization is similar to B u — > 



. The highlighted region is ex- 

TV T . 



The modes D s — > //z/^ and -D s — > rz/ r were shown in [35] to give constraints on (m H +, tan/?) 
in 2HDM II competitive to those from B meson decays. Here we would like to stress the 
additional role played by these observables as tests of the couplings to the second generation of 
quarks in connection with both second and third generation leptons. At tree-level, the D s — > £i/g 
branching ratio is given by [51,52] 



BR{D S 
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with the 2HDM factor 
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(30) 



(31) 



From lattice QCD we adopt the value fn s = 241 ± 3 MeV [50] which has the smallest quoted 
uncertainty. A major source of concern in Eq. fl3Tl) is the dependence on the light quark masses. 
Effects of parametric uncertainties in the masses on the constraints from D s was investigated 
in [35]. We will not repeat this discussion here, but treat the center values as exact quantities. 

Constraints on the 2HDM couplings from the two modes D s — > tv t and D s — ► are 
summarized in Figs. [8] and [U respectively. The excluded regions appear quite similar in these 
figures, as a result of the helicity suppression in the SM, which gives the same relative 2HDM 
contribution to both decay modes. This is interesting for the prospects of testing the hypothesis 
\ TT = A w . A similar test could be performed using the third generation quarks by considering 
B u — > \Wjx in addition to B u — > tv t . Using the same input as above, the SM prediction becomes 
BR(5 U — > fiUfj,) = 4.3 x 10~ 7 . Unfortunately, experiments are not yet sensitive at this level, 
with the current best limit BR(B U — > jju^) < 1.6 x 10~ 6 [53] reported at 95% CL. 

For D s — > one may also combine the result with that from K — > on the coupling 



A, 



iA^, for which we found —40 < \ s 



A w < 30 when m H + = 150 GeV. For m H+ = 300 GeV 
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Figure 8: Constraints on (A CC A TT , X SS X TT ) from BR(D S — > tu t ) for fixed m H + = 150 GeV (left) 
and m H + = 300 GeV (right). The highlighted region is excluded at 95% C.L. The axes are 
normalized similarly to Fig. [61 




Figure 9: Constraints on (A CC A W , A^A^) from BR(D S — * nu^) for fixed m H + = 150 GeV (left) 
and m>H+ = 300 GeV (right). The highlighted region is excluded at 95% C.L. The axes are 
normalized similarly to Fig. [61 
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we obtained essentially no constraint from K — ► jxu^. 

5 Models with a Yukawa Z<i symmetry 

Having discussed constraints on the general model in the previous section, we now specialize to 
the 2HDM with a Yukawa Z2 symmetry. In the language of the general model, this corresponds 
to universal couplings over the fermion generations and types with tan j3 the only free parameter, 
e.g. l/Xu = 1/Acc = —Xbb = —X ss = tan/3 for type II couplings. The coupling patterns in the 
four different models we consider are described in Tabled! The results, as shown in Fig.[TD]in the 
(m#+,tan/3) plane, are superimposed to combine the exclusion regions of all flavor observables 
in the same figure. This also means that some regions of the parameter space are excluded by 
more than one observable, such as the high tan (3 region in the type II model. 

From Fig. [101 we first note the exclusion of low tan (3 < 1 in all four models for ra#+ < 
500 GeV. This exclusion comes as a result of three observables: BK(B — > X s j), A , and AM^. 
The constraints at low tan j3 are similar between the models, since the couplings to the up-type 
quarks are universal. In the type I model, a value of tan/3 > 1 signals decoupling of one 
Higgs doublet from the whole fermion sector. In 2HDM types II and III, which share the same 
coupling pattern for the quarks, there exists a tan /3-independent lower limit of mjj+ > 300 
GeV imposed by BR(I? — » X s j). No generic lower limit on is found in type I and type 
IV models. Constraints for high tan/3 are only obtained in the type II model. The reason 
behind this is that the leptonic and semi-leptonic observables require tan /3-enhanced couplings 
XddXu ~ tan 2 (3 ^> 1 (d — d, s,b) for the contributions to be interesting. In the 2HDM III, and 
IV these couplings are instead always \ddXu. = — 1, while in type I they are proportional to 
cot 2 /3. 

An alternative approach to constraining the parameter space is to assume the 2HDM is 
the correct theory and estimate the parameters by fitting to the available data. To compare 
with our results in Figure [TUl we perform a basic x 2 fit to the flavor data of the parameters 
(ra#+,tan/3) for the Z2 symmetric models. The \ 2 measure is constructed in the usual way 

X 2 = Yli i > where of = (<r° xp ) 2 + (of HDM ) 2 combines experimental and theoretical 

i 

uncertainties in quadrature. Parametric uncertainties are not included. We use the eight 
observables listed in Table [3] as independent probes of the 2HDM. Minimizing the x 2 ■, the first 
observation is that the 2HDM types I, III, and IV, do not contain enough predictive power for 
the flavor observables to restrict the parameter space (m#+,tan/3). In all these models, the 
decoupling limit in one variable (or both) is essentially the best fit. For the 2HDM II on the 
other hand, we obtain a best fit point: m H + = 609 GeV, tan/3 = 5. The fit is illustrated in 
Fig. [TTJ The obtained exclusion region is similar to the combination presented in Fig. [TUJ At 
the level of 2 a, charged Higgs masses m H + < 260 GeV are excluded (m H + < 300 GeV for one 
degree of freedom). The exclusion in tan/3 extends to tan/3 < 1 and tan/3 > 60 for ra#+ = 500 
GeV. Both the lower and the upper limits on tan j3 become stricter for lower rriH+ ■ In Table El 
the results for the 2HDM II are given in more detail. The results for h — > 57 transitions is what 
drives the fit away from decoupling in rriH+- In fact, all other measurements favor new physics 
contributions of opposite sign to those generated by the 2HDM II. This results in a tan (3 close to 
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Figure 10: Excluded regions of the (m H +, tan j3) parameter space for Z 2 -symmetric 2HDM 
types. The color coding is as follows: BR(£> — > X s ^) (red), A _ (black contour), AM# d (cyan), 
B u — > tu t (blue), B — > Dtu t (yellow), K — >• \xv^ (gray contour), D s —> tu t (light green), and 
D s — > fii/p (dark green). 
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Figure 11: Combined parameter estimation in (m H +,tan(3) for the 2HDM II (left). The red 
cross indicates the best fit point ra#+ = 609 GeV, tan/? = 5. The contours show levels of 
Ax 2 = 2.30 (6.18), corresponding to probabilities for la (2 a) Gaussian confidence intervals 
with two degrees of freedom. The center and right panels show variation of A% 2 when changing 
the parameters m#+ and tan (3 around their best fit values. 



its value of minimum influence, tan/5 = y/rrit/mi, ~ 6.5 (tan/3 = \Jm c /m s ~ 3.5) for the third 
(second) generation, since no decoupling limit in tan/3 exists in the 2HDM II. Figure [TT1 shows 
that the \ 2 does not deteriorate significantly when increasing m#+ , illustrating that decoupling 
is not far from being the best fit to the 2HDM II either. The reduced xV d -°- f - = 2 - 1 for the 
2HDM fit, versus x 2 /d.o.f. = 1.8 for the SM reference, confirms that the flavor data does not 
favor a 2HDM II over the SM at this point. 



6 Conclusions 

Allowing for non-universal Yukawa couplings, the 2HDM can modify significantly the theo- 
retical predictions for many flavor physics observables. Using the public codes 2HDMC [27] 
and Superlso [28], we have exploited this fact and used available experimental results to place 
constraints on the couplings of the second and third generation fermions to the charged Higgs 
boson in models with general, flavor- diagonal, Yukawa couplings. Strong and universal con- 
straints are found on the top quark coupling X tt from AM# and b — > sj transitions, requiring 
|Att| < 1 for m>H+ < 500 GeV. Observables such as B u — > tv t or D s — > £ut, with a tree- level 
charged Higgs contribution depending on products of quark and lepton couplings, show sensi- 
tivity at low masses m#+ < 300 GeV both for the second and the third generations. Typically, 
the exclusion limit starts at values [A^A^I 1 / 2 > 10 for m H + = 100 GeV with softer limits for 
higher masses. 

For the versions of the 2HDM with a Yukawa Z2 symmetry, the lower limit m#+ > 300 
GeV from BK(B — > X s '-f) is valid for type II and III models, while no general limit on 
exists for types I and IV. The bound on X tt in the general model translates into a requirement 
of tan/3 > 1. The tree-level observables investigated yield constraints almost exclusively on the 
type II model at high tan/3. For this model, we find that only the observables from b — > S7 
transitions favor a non-zero 2HDM contribution, leading to a best fit point = 609 GeV, 
tan j3 = 5 close to decoupling. 
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Observable 


Experimental 


SM prediction 


AviLr 
AbM 


2HDM fit 


Av 2 

AzrlDM 


Pull 


BR(B - 




3.52 x 10~ 4 


3.07 x 10" 4 


1.65 


3.59 x 10~ 4 


0.04 


0.21 


Anf.B - 

t — nj y 


-> ATM 


3.1 x 10~ 2 


7.8 x 10~ 2 


2.82 


7.0 x 10 -2 


1.88 


1.37 


AMr 


IP / 


0.507 


0.53 


0.08 


0.53 


0.10 


0.32 


BR(B U 


— ► T!/ T ) 

' T J 


1.73 x 10~ 4 


0.95 x 10~ 4 


1.71 


0.95 x 10~ 4 


1.72 


-1.31 


imv 




0.416 


0.30 


0.84 


0.30 


0.84 


-0.91 


Rm{K 


-»■ ^m) 


1.004 


1.000 


0.33 


1.000 


0.33 


-0.58 


BR(D S 




5.8 x 10~ 3 


4.98 x 10" 3 


3.32 


4.98 x 10" 3 


3.36 


-1.83 


BR(D S 


-> rz/ T ) 


5.7 x 10~ 2 


4.82 x 10~ 2 


3.82 


4.82 x 10~ 2 


3.82 


-1.95 








Total x 2 (y)- 


14.6(8) 




12.1 (6) 





Table 3: Best fit of the 2HDM II to the flavor observables compared to the results for the SM. 
The experimental and theoretical uncertainties are given in Table [2j 



The region of 2HDM parameter space where the LHC experiments expect sensitivity (small 
rriH+, large tan/3) [54,55] is very interesting for most of the observables discussed here. Should 
a charged Higgs be discovered at the LHC, the underlying model still needs to be determined. 
A low mass charged Higgs points to either the existence of additional low-energy states (as 
in supersymmetry) , or non-type II couplings for H + . The 2HDM Yukawa sector is accessible 
through flavor physics. In favorable cases, coupling universality between A TT and A MAt can be 
tested using leptonic B and D s decays which are helicity suppressed in the SM. Likewise it 
might be possible to test universality over the fermion generations by combining results from B 
decays with those from K — > \w or D s — > iv. To summarize, a charged Higgs boson discovery 
would provide great opportunities for a rich interplay between flavor and collider physics for 
many years to come. 
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